Phase transition in ultrathin magnetic films with long range interactions 



Laboratorio de Simulagao 



M. RapiniB R.A. DiasQ and B.V. Costfi 
Departamento de Fisica - ICEX - UFMG 30123-970 Belo Horizonte 



MG, Brazil 



O 
O 

(N 

u , 

Oh: 

< 

00 



Oh. 

r 

o 
o 

O 

Oh: 



> 

\o 
o 

o 
\o 
o 

o 

• I— I 



X 



Ultrathin magnetic films can be modeled as an anisotropic Heisenberg model with long range dipo- 
lar interactions. It is believed that the phase diagram presents three phases: A ordered ferromagnetic 
phase (J), a phase characterized by a change from out-of-plane to in-plane in the magnetization (//), 
and a high temperature paramagnetic phase (III)- It is claimed that the border lines from phase 
/ to III and // to /// are of second order and from / to II is first order. In the present work we 
have performed a very careful Monte Carlo simulation of the model. Our results strongly support 
that the line separating phase // and /// is of the BKT type. 



PACS numbers: 



I. INTRODUCTION 



Since the late 80 's there has being; an increasing inter- 
est in ultrathin magnetic films 1, 2, 3, 4, 5, 6]. This inter- 
est is mainly associated to the development of magnetic- 
non-magnetic multilayers for the purpose of giant mag- 
netoresistence applications Q- In addition, experiments 
on epitaxial magnetic layers have shown that a huge va- 
riety of complex structures can develop in the system 
Hi M> H^ 113' Rich magnetic domain structures like 
stripes, chevrons, labyrinths and bubbles associated to 
the competition between dipolar long range interactions 
and strong anisotropics perpendicular to the plane of the 
film were observed experimentally. A lot of theoretical 
work has been done on the morphology and stability of 
these magnetic structures [12, 113, H^- Beside that, it 
has been observed the existence of a switching transi- 
tion from perpendicular to in-plane ordering at low but 
finite temperature [T^ [1^, 113, UM ■ ^^ ^ow temperature 
the film magnetization is perpendicular to the film sur- 
face, as temperature rises the magnetization flips to an 
in-plane configuration. Eventually the out-of-planc and 
the in plane magnetization become zero 19] . 

The general Hamiltonian describing a prototype for a 
ultrathin magnetic film assumed to lay in the xy plane is 

m 



H 



J 2_^ Si ■ Sj — A 2_^ Si 



(1) 



<ij> 



^E 

<ij> 



^i ' ^j _ o 



Si ■ rij j ■ iSj ■ r.ij 



Here J is an exchange interaction which is assumed to be 
nonzero only for nearest-neighbor interaction, D is the 
dipolar coupling parameter, A is a single-ion anisotropy 
and ri,, = r, — fv where r^ stands for lattice vectors. The 
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structures developed in the system depend on the sam- 
ple geometry and size. Several situations are discussed 
in reference I14t| and citations therein. 
Although the structures developed in the system are well 
known the phase diagram of the model is still being stud- 
ied. There are several possibilities since we can combine 
the parameters in many ways. We want to analyze the 
case J > in some interesting situations. A more de- 
tailed analysis covering the entire space of parameters is 
under consideration. 

• Case Z? = 

For £) = we recover the two dimensional (2d) 
anisotropic Heisenberg model. The isotropic case, 
^ = 0, is known to present no transition |2ll| . 
For A > Q the model is in the Ising universality 
class J23| undergoing a order-disorder phase transi- 
tion. If A < the model is in the XY universality 
class. In this case it is known to have a Berezinskii- 
Kosterlitz-Thouless [BKT) phase transition which 
is characterized by a vortex-anti- vortex unbinding, 
with no true long range order Q2t i2^ iM, i2^ ■ 

• Case D^Q 

In this case, there is a competition between the 
dipolar and the anisotropic terms. If D is small 
compared to A we can expect the system to have an 
Ising behavior. If D is not too small we can expect 
a transition of the spins from out-of-plane to in- 
plane configuration [l7|. For large enough D out- 
of-plane configurations become instable such that, 
the system lowers its energy by turning the spins 
into an in-plane anti- ferromagnetic arrangement. 

Earlier works on this model which discuss the phase dia- 
gram were mostly done using renormalization group ap- 
proach and numerical Monte Carlo simulation |l7l llfll 
l2a | . They agree between themselves in the main features. 
The phase diagram for fixed A and J = 1 is schemati- 
cally shown in figure Q in the space {D,T). From Monte 
Carlo (MC) results it is found that there are three re- 
gions labelled in the figure ^ as /, // and ///. Phase / 
correspond to an out-of-planc magnetization, phase // 
has in-plane magnetization and phase /// is paramag- 
netic. The border line between phase / to phase // is 




FIG. 1: A sketch of the phase diagram for the model, (eq. 0. 
Phase I correspond to an out-of-plane magnetization, phase 
II has in-plane magnetization and phase III is paramagnetic. 
The border line between phase / to phase II is believed to be 
of first order and from region I and // to III are both second 
order. 



believed to be of first order and from region / and // to 
/// are both second order. 

Although, the different results agree between them- 
selves about the character of the different regions, much 
care has to be taken because they were obtained by using 
a cut-off, Tc, in the dipolar term. The long range charac- 
ter of the potential is lost, consequently we can expect a 
line of BKT transition coming from region // to region 
///. It is characterized by having no true long range or- 
der. This lack of long range order is prevented by the 
Mermin- Wagner theorem i21f|. The BKT phase transi- 
tion is an unusual magnetic phase transition character- 
ized by the unbinding of pairs of top olog ical excitations 
named vortex-anti- vortex [Ulillillilllllillill. A 
vortex (Anti-vortex) is a topological excitation in which 
spins on a closed path around the excitation core pre- 
cess by 2tt {—2tt). Above Tbkt the correlation length 
behaves as ^ Ri exp(6i-i/2)^ with t={T- Tbkt)/Tbkt 
and ^ — > oo below Tbkt- 

In this work we use MC simulations to investigate the 
model defined by equation ^ We use a cutoff, Tc, in the 
dipolar interaction. Our results strongly suggest that 
the transition between regions // and /// is in the BKT 
universality class, instead of second order, as found in 
earlier works. 



spacing, as suggested in the work of Santamaria and co- 
workers [l7| . 

We have performed the simulations for temperatures 
in the range 0.3 < T < 1.2 at intervals of AT = 0.1. 
When necessary this temperature interval is reduced to 
AT ~ 0.01. For every temperature the first 5 x 10^ MC 
steps per spin were used to lead the system to equilib- 
rium. The next 10^ configurations were used to calcu- 
late thermal averages of thermodynamical quantities of 
interest. We have divided these last 10^ configurations 
in 20 bins from which the error bars are estimated from 
the standard deviation of the averages over these twenty 
runs. The single-site anisotropy constant was fixed as A 
= 1.0 while the D parameter was set to 0.10,0.15 and 
0.20. In this work the energy is measured in units of J S^ 
and temperature in units of JS^/ks, where fcs is the 
Boltzmann constant. 

To estimate the transition temperatures we use finite- 
size-scaling (FSS) analysis to the results of our MC sim- 
ulations. In the following we summarize the main FSS 
properties. If the reduced temperature ist = {T — Tc)/T, 
the singular part of the free energy, is given by 



T(L,T) = L-(2-")/^jr(iii/'') 



(2) 



Appropriate differentiation of F yields the various ther- 
modynamic properties. For an order disorder transition 
exactly at Tc the magnetization, M, susceptibility, ^ and 
specific heat, C, behave respectively as |3lll32|. 



(3) 



In addition to these an important quantity is the forth 
order Binder's cumulant 
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(4) 



3 < m"^ >^ 

For large enough L, curves for Ui{T) cross the same 
point U* at T = Tc- For a BKT transition the quan- 
tities defined above behave in a different way. Due to the 
Mermin- Wagner theorem there is no spontaneous mag- 
netization for any finite temperature. The specific heat 
present a peak at a temperature which is slightly higher 
than Tbkt- Beside that, the peak height does not de- 
pends on L. Because models presenting a BKT transi- 
tion have an entire critical region, the curves for Ui{L) 
just coincide inside that region presenting no crosses at 
all. Below we present MC results for three typical re- 
gions. When not indicated the error bars are smaller 
than the symbol sizes. 



II. SIMULATION BACKGROUND 

Our simulations are done in a square lattice of volume 
L X L (L = 10, 20, 30, 40, 50, 80) with periodic boundary 
conditions. We use the Monte-Carlo method with the 
Metropolis algorithm [HlllllSlp. To treat the dipole 
term we use a cut-off Vc = 5a, where a is the lattice 



III. SIMULATION RESULTS 
A. Case D=0.1 

For D = 0.1 we measured the dependence of the out-of- 
plane magnetization, M^, and the in-plane magnetization 



, Mxy, as a function of temperature for several values of 
L (See figure 121). The figures indicate that in the ground 
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FIG. 2: Out-of-plane (left) and in-plane (right) magnetization 
for D = 0.1. The ground state is ferromagnetic. There is no 
in-plane spontaneous magnetization. 

state the system is aligned in the z direction. Approxi- 
mately at T « 0.70 the Mz magnetization goes to zero, 
which gives a rough estimate of the critical temperature. 
The in-plane magnetization has a small peak close to 
T « 0.70. However, the height of the peak diminishes as 
L grows, in a clear indicative that it is a finite size arti- 
fice. The behavior of the specific heat, susceptibility and 
Binder's cumulant, are shown in figures|31E]and|31respec- 
tively. The results indicate a order-disorder phase tran- 
sition in clear agreement with references [la, llfl lOi H^ ■ 
The vortex density in the xy plane (Figure ^ has a 
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FIG. 3: Specific heat as function of temperature for D — 0.1. 

very shallow minimum near the estimated critical tem- 
perature which is almost independent of the lattice size. 
The growth of the number of vortices when the tempera- 
ture is decreased is related to the disordering in the plane 
when the magnetic moments tend to be in the z direc- 
tion. We have performed a finite size scaling analysis 
of the data above by plotting the temperature T^ as a 
function of 1/L for the specific heat, the susceptibility 
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FIG. 4: Out-of-plane susceptibility as function of temperature 
for D = 1 





1 1 1 1 1 1 1 1 1 1 1 1 1 1 


' 1 ' 1 ' 1 ' 


0.6 


W* 


»-oL=10 
»-oL=20 




|U 


»-«L=30 


05 


i-iL=40 




iflT k 


<m=50 


0.4 




^^L=80 


0.3 


ffii 
I'iii 

Til 


1 


0.2 


i 


* 


0.1 


>i\^ 


*- 




wV- 


•v '»--» . 





1 1 1 1 1 1 1 1 1 1 1 1 1 1« 


T 1 1 1 1 1 1 1 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 
T 

FIG. 5: Binder's cumulant as function of temperature for 
D = 0.1. 
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FIG. 6; Vortex density in the xy plane for D — 0.1. 



and the crosses of the forth-order cumulant. The results 
are shown in the table QJ By linear regression we have 
obtained the critical temperature as T^°° = 0.682(2). An 
analysis of the behavior of the maxima of the specific 
heat, Craax^ (See figure [T7|) as a function of the lattice 
size shows that it behaves as Cmax oz InL, indicating a 
second order phase transition. In the phase diagram we 
crossed the second order line labelled c. 
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TABLE I: Critical temperature T^ of the specific heat, C, 
susceptibility, x, and the crosses of the fourth-order Binder's 
cumulante Ui as a function of the lattice size L. Data are for 
D = 0.10 



Ti « 0.35. The second peak appears at T2 « 0.65 and 
seems to be independent of the lattice size. In the 
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FIG. 8: Specific heat for D = 0.15. 



B. 



Case D=0.15 



In this region of the parameters, it was observed a 
transition from an out-of-plane ordering at low temper- 
atures to an in-plane configuration as described by the 
magnetization behavior shown in Fig. [T] We show Mz 
and Mxy in the same figure for comparison. The out-of- 
plane magnetization goes to zero at T « 0.35 while an 
in-plane magnetization sets in. This phenomenon has al- 
ready been reported experimentally [iSl [Ig and it is due 
to the competition between the easy axis anisotropy and 
the dipolar interaction. The specific heat curve presents 
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FIG. 7: Mz and M^y (open and full symbols respectively) for 
L> = 0.15. 

two peaks (See figure|Hl). The peak at low temperature is 
pronounced and is centered in the temperature in which 
occurs the rapid decrease of the in-plane magnetization. 



FIG. 9: In-plane (left) and out-of-plane (right) susceptibility 
for D = 0.15. 

figure ini we show the in-plane and out-of-plane suscepti- 
bilities. The out-of-plane susceptibility presents a single 
peak close to Ti « 0.35. The in-plane susceptibility has 
a maxima at T2 « 0.65 beside the peak at Ti, indicat- 
ing two phase transitions. The Binder's cumulant for 
the in-plane and out-of-plane magnetization are shown 
in figures ^1 Except for the case L = 10 the curves for 
different values of the lattice size do not cross each other 
indicating a BKT transition at T « r2- Beside that, the 
in-plane cumulant has a minimum at T « Ti, which is a 
characteristic of a first order phase transition 31, ^Z|. 

The vortex density is shown in figure [TTl Its behavior 
is similar to that one shown in figure In the phase 
diagram we crossed the region I to II (Ti) and II to III 
(T2). The maxima of the specific heat are shown in figure 
1171 as a function of L. It is clear that after a transient 
behavior it remains constant indicating a BKT transi- 
tion. A FSS analysis of the susceptibility (see table |n] 
) gives the BKT temperature T'g'i^j, = 0.613(5). In the 
phase diagram we crossed the first order line labelled a 
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FIG. 10; In-plane (left) and out-of-plane (right) Binder's cu- 
mulant as function of temperature for D = 0.15. Observe that 
the in-plane cumulant has a minimum at T ~ 0.35 indicating 
a first order phase transition. After the minimum the curves 
do not cross each other having the same behavior (Except the 
spurious case L = 10.) up to T ~ 0.65 when they go apart. 
That is an indication of a BKT phase transition. 
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FIG. 11: Vortex density as function of temperature for D 
0.15. 



(Ti) and the line labelled b {T2). 



C. Case D=0.20 

In figure lT^ we show the in-plane and out-of-plane mag- 
netization curves for several lattice sizes and D = 0.20. 
We observe that as the lattice size L goes from L = 10 to 



L = 80, both magnetization decrease. It can be inferred 
that as the system approaches the thermodynamic limit, 
the net magnetization should be zero. Therefore, the sys- 
tem does not present finite magnetization for any tem- 
perature T ^ 0. The specific heat (Figure [T^ presents 
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FIG. 12: Mxy and Mz (open and full symbols respectively) 
for D = 0.2. 

a maximum at T w 0.75. The curves are for different 
values of L. We observe that the position of the maxima 
and their heights are not affected by the lattice size, all 
points falling in the same curve. 
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FIG. 13: Specific heat for D = 0.2. The line is a guide to the 
eyes. 



T/- 0.729 0.69^ 
L 10 20 



0.678 0.670 0.650 0.638 
30 40 50 80 



TABLE II: Critical temperature T^ as a function of the linear 
size L for the susceptibility x. 



In the figure ini we show the in-plane and out-of-plane 
susceptibilities respectively, x^^ does not present any 
critical behavior, x^^ presents a peak which increases 
with L. For the Binder's cumulant there is no unique 
cross of the curves. (Except for the L = 10 curve, which 
is considered too small to be taken in to account.). This 
behavior indicates a BKT transition at Tbkt ~ 0.63. 
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FIG. 14: In-plane and out-of-plane susceptibility for D — 0.2. 
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The vortex density, shown in figure [T^ is almost indepen- 
dent on the lattice size. In addition, we did a FSS analy- 
sis of the susceptibility (see table IIIl|) and the maxima of 
the specific heat. The specific heat is shown in figure [T7I 
Its behavior indicates a BKT transition. The analysis of 
the susceptibility gives Tgj^j. = 0.709(5). In the phase 
diagram we crossed the line labelled b. In our results 
we could not detect any other transition for D = 0.20, 
indicating that: The line labelled a ends somewhere in 
between 0.15 < D < 0.20 or the crossing at a occurs at 
a lower temperature (T < 0.30) outside the range of our 
simulated data. 



T^ 0.829 0.781 0.768 0.753 0.750 0.729 
L 10 20 30 40 50 80 



TABLE III: Critical temperature T^ as a function of the lin- 
ear size L for the susceptibility x- 
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FIG. 16: Vortex density in the xy plane for D = 0.2. 



5 


' 1 


1 


1 ' 1 


' 1 


4.5 


- 






^ 


4 


- 




0^^'^ 


- 


3.5 


- 






- 


C 3 


^^ 






- 


2.5 


- 






- 


? 


_ 






_ 






D 






1.5 
1 




1 




1 


1 , 1 


1 




. 2.5 


3 


3.5 4 
InL 


4.5 5 



FIG. 17: Maxima of the specific heat as a function of the 
lattice size. 



IV. CONCLUSIONS 

In earlier studies several authors have claimed that the 
model for ultrathin magnetic films defined by the equa- 
tion n presents three phases. Referring to figure Q] it is 
believed that the line labelled a is of first order. The line 
b and c are of second order. Those results were obtained 
by introducing a cut off in the long range interaction of 
the hamiltonian. In the present work we have used a 
numerical Monte Carlo approach to study the phase di- 
agram of the model for ,/ = A = 1 and D ~ 0.10,0.15 
and 0.20. In order to compare our results to those dis- 
cussed above we have introduced a cut off in the long 
range dipolar interaction. A finite size scaling analysis 
of the magnetization, specific heat, susceptibilities and 
Binder's cumulant clearly indicates that the line labelled 
a is of first order and the line c is of second order in 
agreement with other results. However, the b line is of 



BKT type. After analysing the results obtained, some 
questions come out: 

1. Is it possible the existence of a limiting range of 
interaction in the dipolar term beyond which the 
character of the transition changes from BKT to 
second order ? 

2. How does the line labelled a end in the phase dia- 
gram ? 

3. What is the character of the intersection point of 
the three lines in the phase diagram ? 

In a very preliminary calculation Rapini et al.|l8J stud- 
ied the model with true dipolar long range interactions. 
Their results led them to suspect of a phase transition of 



the BKT involving the unbinding of vortices-anti- vortices 
pairs in the model. However, to respond those questions 
it is necessary to make an more detailed study of the 
model for several values of the cut off range re- In a sim- 
ulation program we have to be careful in taking larger Tc 
values since we have to augment the lattice size propor- 
tionally to prevent misinterpretations. 
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